Coupled Mode Effects on Energy Transfer Rates in Two Temperature Plasmas 
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We investigate the effects of collective modes on the temperature relaxation rates in fully coupled 
electron-ion systems. Firstly, the well-understood limit of weakly coupled plasmas is considered 
and the coupled mode formula within the random phase approximation is derived starting from the 
Lenard-Balescu kinetic equation. We show how the frequency integration can be performed by stan- 
dard methods without applying additional approximations. Due to the well-defined approximation 
scheme, the results can serve as a benchmark for more approximate theories and numerical simu- 
lations in this limit. The coupled mode electron-ion transfer rates show a considerable reduction 
compared to the Fermi-Golden Rule approach for certain parameters and very small changes for 
other systems. We demonstrate how these coupled mode effects are connected to the occurance of 
ion acoustic modes and under which conditions they occur. Interestingly, coupled mode effects can 
also occur for plasmas with very high electron temperatures; a regime, where the Landau-Spitzer 
approach is believed to give accurate results. Finally, we extend the approach to systems with 
strongly coupled ions by applying static local field corrections. This extension can substantially 
increase the coupled mode effects. 

PACS numbers: 52.25.Dg, 52.25.Kn, 52.27.Gr 
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I. INTRODUCTION 

Experimental characterization of warm dense matter 
and strongly coupled plasmas becomes increasingly pow- 
erful with the development of new methods for the cre- 
ation and probing of such states. For instance, x-ray scat- 
tering allows not only for the measurement of equation 
of state data (like in traditional shock wave experiments, 
see, e.g., Refs. [H, HI), but also to obtain structural, dy- 
namic, and collective properties of matter [3, EL IE IE 0] • 
With these new possibilities, one is now able to probe 
the physics of high energy density matter as it is encoun- 
tered during inertial confinement fusion or in the interior 
of planets. 

The creation of the these high energy density states in 
the laboratory requires large and fast energy inputs into 
matter. Since static techniques, like diamond anvil cells, 
are restricted to lower temperatures and densities by the 
given stress point of the material, one relies on dynamic 
experiments applying intense particle beams and lasers to 
heat and compress the material under investigation. In- 
evitably, highly nonequilibrium states are produced with 
the energy being pumped mainly either into the ion or 
the electron subsystem. Such states are of great interest 
by itself since they reveil many dynamic and collective 
properties of the system more clearly. 

A good understanding of nonequilibrium states is also 
needed for definitive measurements of equilibrium prop- 
erties. After a very short time of approximately an in- 
verse plasma frequency electrons and ions have estab- 
lished temperatures within their subsystems and a full 
kinetic description is not necessary later on. Apart from 
the final hydrodynamic response, temperature equilibra- 
tion thus takes the longest time of all relaxation pro- 
cesses driven by the intial energy deposition and defines 



the minimum time delay between the pump and probe 
pulses needed for equilibrium measurements. 

Temperature equilibration is furthermore interesting 
since its time scale is experimentally accessible. Indeed, 
experiments investigating dense plasmas, both laser- 
and shock-produced, found relaxation times considerably 
longer than predicted by classical Landau-Spitzer for- 
mula [1, Q. Two shortcomings of this easy-to-use ap- 
proach were associated with these discrepancies: the ne- 
glect of collective excitations and the use of classical col- 
lisions. A full quantum binary collision approach yields 
however even larger energy transfer rates [lol [Til . [r3 |. 
Considering independent collective modes in the elec- 
tronic and ionic subsystems within the Fermi-Golden- 
Rule (FGR) approach yields rates very close to the 
Landau-Spitzer results [l3(. Therefore, collective modes 
in fully coupled electron-ion systems seem to be the only 
candidate to explain the lower electron-ion energy trans- 
fer measured. 

More complex descriptions of temperature relaxation 
involve the interplay of all contributions to the internal 
energy. In addition to the kinetic parts interesting for 
temperature relaxation, correlations and exchange ener- 
gies [3 HE 16 1 and the ionisation equilibrium including 
excitations 17, [IE Ell can be important. The latter 
processes have been shown to considerable influence the 
relaxation under conditions far from equilibrium, but are 
not sufficient to explain the measured slower temperature 
equilibration. 

Here, we assume ionisation, recombination, and charge 
transfer processes to be completed and concentrate on 
the electron-ion energy transfer including collisions and 
collective excitations in a fully coupled system. We use a 
quantum statistical description which avoids any ad hoc 
cutoffs known from classical descriptions and enables us 
to rigorously derive formulas for the energy transfer rates 
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within a given approximation scheme. 

The occurance and magnitude of coupled mode effects 
on temperature relaxation are still under discussion (see 

Refs. HUHHHHHlii]). In the be S innin S' we therefore 
consider weakly coupled plasmas, where the classical cou- 
pling parameter, Tj = Z 2 e 2 joj ksTj with Oj — (3/47rri ? ) 1 / 3 
is small for all species. In this limit, we can employ the 
well-established random phase approximation (RPA) to 
describe the dynamic response of fully coupled systems 
and the quantum version of the Lenard-Balescu equation 
[25l l26l . [27j on the kinetic level. We show that the weak 
coupling version of the coupled mode formula derived by 
Dharma-wardana & Perrot [2(| follows. Moreover, we re- 
write the density response functions of their approach in 
terms of dielectric functions and obtain a form that can 
be evaluated by standard integration procedures. 

Under certain conditions, we observe a lowering of the 
energy transfer rates compared to the Landau- Spitzer 
and the FGR approaches in the order of a factor of two. 
These differences are connected to the occurance of ion 
acoustic modes and a related redistribution of weight in 
the dynamic density response toward smaller frequencies. 
This fact also explains why a simpler approximation like 
FGR is suitable for other parameters. The most striking 
fact is that coupled mode effects can occur also for very 
high electron temperatures where the Landau-Spitzer ap- 
proach is believed to work well. 

We have good control over our approximation scheme. 
The results in the well-understood weak coupling limit 
can thus serve as a benchmark for more approximate ap- 
proaches. In particular, numerical simulations such as 
applied in Refs. [H, [2^] should be tested against these 
analytic results since they are based on classical mechan- 
ics and apply pseudo-potentials to approximately incor- 
porate the qantum nature of the electrons as an uncon- 
trolled approximation. 

As an extension to the RPA scheme, we consider 
strongly coupled ions where the additional correlation 
effects are treated on the level of static local field correc- 
tions. The additional shifts in the mode structure can 
further reduce the electron- ion energy transfer. For com- 
pressed silicon as discussed in Ref. @, we report large 
reductions although coupled mode effects on the RPA 
level do not change the rates much. In addition, we dis- 
cuss the influence of degeneracy, mass and ion charge. 



II. ENERGY TRANSFER RATES 

A. General Kinetic Description 

Changes of the kinetic energy of species a, that is E a , 
are fully determind by the changes of the one-particle 
Wigner distribution of the species f a 
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In the second line, two-particle collision integrals l ab were 
introduced by using a general kinetic equation for homo- 
geneous and isotropic systems: df a /dt = ^ b I a b- The 
equations for the electron and ion species are coupled by 
the collision integrals I e i — Ii e and, thus, the total energy 
of the plasma is conserved. 

There exists a hierachy of approximations for the col- 
lision integrals I a b [13] which via Eq. (fTJ) also defines the 
quality of the related energy transfer rates. The differ- 
ent kinetic equations can be devided into two classes: 
the first describes the interactions by binary collisions 
only (Landau equation (30j and Boltzmann-like equa- 
tion HH), the second includes also collective effects such 
as dynamic screening (Lenard-Balescu-type kinetic equa- 
tions [ID, [26[ ) . The mutual influence of electrons and ions 
via the common dielectric function of the plasma (cou- 
pled collective modes) is here included and predicted to 
strongly modify the energy transfer in two-temperature 
plasmas |20l [21j|. 



B. Classical and Quantum Binary Collisions 

The easiest approximation for the two-particle colli- 
sion integral considers classical binary collisions. The 
corresponding energy transfer rates were derived in the 
first description of temperature relaxation by Landau and 
Spitzer Hfl H2 



8 F ls 



3 , Ti — T e 
—n e KB 

2 Tf>7* 



with the clcctron-ion relaxation time 
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The plasma properties define the Coulomb logarithm 
In A, where A is the ratio of the maximum and mini- 
mum impact parameters. Here, we use a Coulomb loga- 
rithm of the form lnA = 0.51n(I + \ 2 e /{o]_ + \ 2 dB )) with 
the electron screening length A e = (fcsT e /47re 2 n e ) 1 / 2 , 
the distance of closest approach g± = Z b e 2 /m e v 2 h , the 
deBroglie wave length \<ib = Ji/meVthi an d the thermal 
velocity vth = (kBTe/m,,) 1 / 2 . This form, which follows 
by considering hyperbolic orbits of the electrons, has the 
advantage to give non-negative results even for the dense 
plasmas considered in this paper. For a more extensive 
discussion see Ref. (loj . 

The Landau-Spitzer formula ([3]) suffers from a number 
of deficits. The most crucial one is the classical descrip- 
tion of clcctron-ion collisions which may result in neg- 
ative Coulomb logarithms and rates. Well-defined en- 
ergy transfer rates in the binary collision approximation 
can however be derived based on the collision integral of 
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the quantum Boltzmann equation [lfj, [34{ . If the colli- 
sion cross sections needed are calculated from the two- 
particle Schrodinger equation, these rates consider strong 
electron-ion collisions as well but collective excitations 
cannot be included since statically screened interactions 
are used. Most binary collision approaches are also lim- 
ited to systems with nondegenerate electrons due to the 
way the collision cross sections are calculated. 



C. Collective Modes in Weakly Coupled Plasmas 

The collective excitations of a weakly coupled plasma 
are fully accounted for in the well-known random phase 
approximation (RPA). Since the electrons and ions are 
treated as one combined system in this approach, their 
mutual influence is naturally accounted for and, so called, 
coupled collective modes arise. This becomes particularly 
clear when considering the dielectric function on the RPA 
level that includes a sum over all species 



.RPA 
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where xla i s ^ ne f pee density response function for the 
particles of species a. 

The kinetic equation corresponding to the RPA is the 
Lenard-Balescu equation [25|, [26J. Here, we will use 
its quantum generalization that accounts for degeneracy 
[27l | . The energy transfer rates are again obtained by in- 
serting the appropriate collision integral into the energy 
balance |T]). Details of the derivation can be found in 
Appendix fA] and the final result is 
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(5) 

n B (u) — [exp(huj/kBT a ) — is a Bosc function that 
characterises the occupation number of the collective 
modes and V e i(k) = 4:irZe 2 /k 2 is the pure Coulomb po- 
tential. The zeros of the total dielectric function in the 
denominator determine the collective excitations of the 
fully coupled system. Thus, the energy transfer rate ([5]) 
carries the label 'CM' for coupled mode. It is applica- 
ble for weakly coupled plasmas without any restriction 
with respect to degeneracy. For weakly coupled plas- 
mas, equation ([5]) is equivalent to the CM formula given 
by Dharma-wardana & Perrot (see Eq. (50) in Ref. (20j ) 
which is demonstrated in Appendix 151 

We would like to emphasis again that expression ([5]) 
represents a CM formula for the energy transfer rates 
in weakly coupled, two-temperature systems. It was 
derived from the well-established (quantum) Lenard- 
Balescu equation in RPA without further approxima- 
tions. Therefore, it can serve as a benchmark for more 



approximate approaches and numerical simulations in 
the weak coupling limit. 

One of the more approximate, but easier to calculate, 
approaches for the energy transfer rates is the Fermi 
golden rule (FGR) approach. Here, electrons and ions 
are treated as separate subsystems. Accordingly, both 
subsystems have their own distinct dielectric functions 
e a PA = VaaXaa- The expression for the energy trans- 
fer rates reflects this fact by having a product of these 
dielectric functions in the denominator 120, [33| 
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It can be shown that the FGR rate reduces to Landau- 
Spitzer (LS) if classical expressions for the dielectric func- 
tion are used [13l . [33| . 

Although both, the CM expression (© and the FGR 
formula ([6]), contain effects of collective excitations, they 
describe quite different systems since the spectra of the 
modes do not agree. To evaluate the corresponding dif- 
ferences in the energy transfer rates is the main purpose 
of this paper. 



D. Incorporation of Strong Coupling Effects 

In dense plasmas, the Coulomb correlations between 
the particles can strongly effect the mode spectrum and 
scattering cross sections. The extensions towards strong 
coupling in the CM formula were already given in the 
original derivation by Dharma-wardana & Perrot [2p| . 

Strong electron-ion collisions may be accounted for by 
a weaker electron- ion pseudo-potential: V e j(fc) — ► U e i(fc). 
Spatial correlations, in particular in the ionic subsystem, 
can be described by a dielectric function that includes 
local field corrections. Here, we will use static local field 
correction (SLFC) for the strongly coupled ions only. In 
this case, the ionic dielectric function reads [35| 
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The local field factor Gu(k) is connected to the ionic 
structure factor Sa(k) via 



G u {k) = 1 



knTi 



niVu{k) \Su(k) 



1 



(8) 



The spatial correlations in Su(k) can be calculated by 
molecular dynamics or Monte Carlo simulations as well 
as integral equations like the hypernetted chain approach 

It should be mentioned here that the changes due to 
SLFC influence only the CM rates. Within the FGR ap- 
proach, the w-integration is essentially determined by the 
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FIG. 1: (Color online) Dispersion relations, Ree(fc, uo) = 0, for 
the modes in a hydrogen plasma with an ion temperature of 
Ti — 10 4 K, a density of n = 10 22 cm -3 , and different elec- 
tron temperatures. Frequencies are normalized to the elec- 
tron plasma frequency. In the FGR approach, there exists an 
electron and an ion plasmon mode (upper blue dotted line) 
whereas all ion modes become acoustic in the coupled system. 



f-sum rule [13, [33] which gives the same result indepen- 
dent of the approximation level used. Therefore, strong 
coupling effects are always included in the FGR rates. 



III. EVALUATION OF THE COUPLED MODE 
AND FGR EXPRESSIONS 

The occurance of collective modes in the systems is 
connected with zeros in the dielectric function and results 
in sharp peaks in the ^-integration. These peaks prohibit 
the evaluation of Eqs. §5$) and ([6]) by standard integration 
routines. We will first show under which conditions these 
modes occur and then describe the numerical integration 
procedure. 



Mode Structure in CM and FGR 



The ionic modes used in the FGR and CM approaches 
are qualitatively different. For small momenta, the ion 
subsystem in FGR has always a plasmon mode which 
is similar to the electron plasmon (just scaled by Zi 
and mi) as well as an acoustic mode. The coupling be- 
tween the electronic and ionic subsystems as included in 
the CM formula leads to screening of the ionic plasmon 
mode which therefore turns into a weakly damped up- 
per acoustic branch of the ionic spectrum. The finite 
mode frequency at k = is thus shifted to much smaller 
frequencies that linearly increase with k. Interestingly, 
this new acoustic branch becomes more plasmon-like for 
larger temperature differences and momenta and coin- 
cides with the prediction for an isolated ion system in the 
unscreened limit with T e — >oo. For the important small 
k values, the modes are however for all finite electron 
temperatures acoustic (in contrast to the FGR modes). 

Fig-fflalso clearly demonstrates that the ionic modes in 
the coupled system cease to exist if the temperature dif- 
ference, T e —Ti, becomes too small. The mode frequency 
is, for small k, in very good approximation given by the 
zeros of the real part of the dielectric function. While 
the product of dielectric functions in the FGR approach 
has always zeros, the summation over species in the full 
dielectric function J3J puts strict limits to the occurrence 
of ionic modes and changes their location in w-space; a 
fact that is crucial for the understanding of the CM effect 
on the energy relaxation. 

To estimate under which conditions ion acoustic modes 
in the coupled plasma occur, we have to balance the 
electronic terms which can be treated statically and the 
frequency-dependent ionic part. It is well-known that a 
zero in the real part of the full dielectric function can be 
achieved under the condition Ti <C T e [37| . A more pre- 
cise analysis (see Appendix [C]) shows that either of the 
relations [39| 



T t < 0.27 ■ Zi T e 
Ti < 0.27 ■ Zi Tp 



for n e Ag < 1 
for n e Ag > 1 



(9) 



must hold to allow for an ion acoustic mode in coupled 
electron-ion systems. The upper case gives the limit for 
nondegenerate electrons while the second line sets the 
limit for highly degenerate electrons (A 2 = 2nh/m e kBT e 
is the thermal wave length). In the high density case, 
the important temperature scale is thus defined by the 
electron Fermi temperature Tp — ?i 2 (37r 2 n e ) 2 / 3 /2/esm e . 
Accordingly, sharp ion acoustic modes may also occur 
for T e <T if the electron density is high enough. 



The electron modes have the well-known acoustic and 
plasmon-like branches and are almost unaffected by the 
presence of the ions. While the acoustic branch is 
strongly damped, the plasmon results in a sharp peak at 
the plasmon frequency oj e ~ (<liTe 2 n e /m e ) 1 ' 2 (see Fig.QJ. 
However, these modes are practically not included in the 
^-integration which is effectively limited by the function 
ImXii{k,uj) to frequencies less than the ion plasma fre- 
quency uoi « (47rZ 2 e 2 n i /mj) 1 / 2 . 



B. Evaluating the FGR Formula 

The FGR formula for the energy transfer rates (O al- 
ways features sharp peaks at small wave numbers that are 
roughly located at the ion plasmon frequency. The exact 
value of the mode frequency f2 can be easily obtained 
numerically. To avoid the sharp peaks, these contribu- 
tions around the location of the mode are cut out. The 
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cut out contributions can be easily calculated from the 
differences to an w-integral that allows for an evaluation 
by the f-sum rule j4(|. Thus, we have to compute 



du> Ime e (w) Im£i(w) 
— OJ ATVei (UJ ) - -o- -n- 

2n keH k,M 



0.001 



du> Im£j(u;) 
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Here, AN e i(uj) — n e B (uo) — n B (uo) is the difference in the 
mode occupation numbers. The function F(u>) includes 
the electronic dielectric functions and AN ei (to), that is, 
F(u) = A7V ei (w) Ime e (aj)/|e e (ci;)| 2 . In this way, the sharp 
plasmon peak has been removed and then treated exactly. 
All remaining integrals are over smooth functions. 



C. Evaluating the CM Expression 

Unfortunately, an evaluation as for the FGR formula 
is impossible for the CM integrand. Here, the electronic 
and ionic parts do not separate and one cannot expand 
the integrand into a full dielectric function needed for the 
f-sum rule. 

The functions under the w-integral can, however, be 
rearranged to allow for an integration with standard rou- 
tines. For this aim, we artificially decompose the total 
dielectric function into an electron and an ion part 

e RPA {k,u) = l+[Ree e (k,uj) - 1] + Ime e (fc,w) 

+ [Reei(k, u) - 1] +Ime;(fe,u;) . (11) 

Then, we express the density respose functions in terms of 
dielectric functions: Im e RPA (k,w) = V aa (k) Im x° a (k,u)). 
In the case of a Coulombic interaction, Vj = VuV ee holds 
and can be used to do the transformation above. The uj- 
integral has now the form 
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This form can still have very sharp peaks at the posi- 
tions of the weakly damped upper ion acoustic modes. 
However, these peaks are limited in height since the real 
part of the total dielectric functions vanish at the mode 
frequency and we have 



lim 

Ree^O 



Ime e (fc,Lij) Im£i(fc,w) 



-RPA 



(*, W )r 



Im£ e (fc,u) Im£i(fc,u;) 
|Ime e (fc, ui) + Im£i(fc, u)\' 



<1. (13) 



This rearrangement makes a brute force approach for in- 
tegrating the CM equation (fT2|) feasible since extremely 
narrow peaks have now no contribution to the integral. 
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FIG. 2: (Color online) Examples the combination of dielectric 
functions in the ui- integrand of the CM and FGR expressions, 
i.e., fcM = Ime e Imei/|ej 2 and /fgh = Im£ e Im£i/|£ e | 2 |£i| 2 , 
respectively. Considered is a small wave number of k — a^ 1 
and a hydrogen plasma with Ti = 10° K, T e = 10 6 K, and 
n — 10 24 cm" 3 . Note: the height of the ion plasmon peak in 
the FGR approach is 1500 a.u. 



D. Results of the ^-integration 

Example of the w-integrands for the two approaches 
are plotted in Fig. [2] The wave number has been cho- 
sen small enough that the ion acoustic mode exists in 
the CM approximation. It is clearly visible that the 
ion mode is shifted to lower frequencies. Moreover, the 
height of this peak has reduced dramatically indicating 
much larger damping. Although the particle excitation 
(left broad peak) is increasing, it cannot compensate the 
loss in weight. For a pure ion system, as in the FGR 
approach, almost the entire contribution stems from the 
plasmon peak whereas the main contribution to the CM 
integral comes from the particle peak. 

The above points are the basics to understand the 
lowering of the energy transfer due to coupled collec- 
tive modes. Figure [3] demonstrates that the CM effect is 
caused by differences at small wave numbers. Here, the 
ion acoustic mode exists and the influence of the elec- 
trons shifts the weight to smaller k as shown in Fig. [T] 
Thus, the overall contributions to the CM integral are of 
the same order as the particle peak of the FGR which 
itself is negligible compared to the ion plasmon peak. 
Accordingly, small fc-values give almost no contribution 
in the CM expression compared to the FGR results as 
shown in Fig. [3l If these small k are important, sig- 
nificantly reduced CM rates, compared to FGR, can be 
expected. For large momentum transfers k, the situation 
is different: here, the ion acoustic mode (in CM) and the 
plasmon mode (in FGR) do not exist anymore and the 
w-integrals of both approaches merge. If such larger k 
dominate in expressions (O and ((6j> , no coupled mode 
effects can be expected. 
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FIG. 3: (Color online) Examples for the w-integral in the CM 
and FGR expressions for hydrogen as in Fig. [2] 



IV. RESULTS AND DISCUSSION OF THE 
ENERGY TRANSFER RATES 



FIG. 4: (Color online) Electron-ion energy transfer rates in 
fully ionized hydrogen with n = 10 26 cm -3 and Ti = 10 5 K 
versus electron temperature. The applied approximations are 
the coupled mode formula within RPA ([5} and with static 
local field corrections (|8]), Fermi's golden rule approach ((6|, 
and the Landau-Spitzer formual ([2]). 



We now turn to the analysis of energy transfer rates 
in dense two-temperature plasmas applying the CM, the 
FGR and the LS approaches. We consider a variety of sit- 
uations typical for laser- as well as shock-produced plas- 
mas. First, we will concentrate on weakly coupled plas- 
mas and describe the plasma within the random phase 
approximation. Based on the methods presented above, 
this allows us to calculate the CM energy transfer rates 
on a well-established approximation level (similar to that 
of the quantum Lenard-Balescu equation) without any 
further approximations. In a second step, we will include 
static local field corrections to extend the CM approach 
to plasmas with strongly coupled ions. 



A. Results for Hydrogen within RPA 

In a first example, hydrogen at a high density, as 
needed for the simulation of inertial confinement fusion 
scenarios, is considered in Fig. [¥] In this example, we 
have relatively warm ions and hotter electrons (the range 
of T e has been extended to see the limiting behavior). All 
approaches show similar curves: i) a strong increase of 
the rates when the electron temperature is just a few 
times the ion temperature, ii) maximum of the energy 
transfer, and iii) a LS-like asymptotic reduction of the 

1/2 

energy transfer ~ T e . However, considerable quantita- 
tive differences arise when comparing the approximation 
levels presented. 

In the first region, the FGR, CM and CM+LFC results 
agree since here ion acoustic or coupled collective modes 
do not existent. As expected, the LS formula gives dif- 
ferent results in this region since the Coulomb logarith 
is not well-defined for these plasmas with relatively cold, 
degenerate electrons. At electron temperatures of about 



T e = 10 6 K, ion accoustic modes start to occur at small k 
and, accordingly, the CM results show increasing devia- 
tions from the FGR rate. 

Around the maximum of the energy transfer rates the 
ion acoustic modes are fully developed. Since they result 
in smaller w-integrals for small k (see Fig. [3]), the CM 
results show here a considerable lowering compared to 
FGR rates. Although strong coupling effects are small 
in this example, the effect of the local field corrections is 
also largest in this region. 

For very high electron temperatures, the Landau- 
Spitzer formula is believed to be accurate due to the very 
weak electron-ion and electron-electron coupling. How- 
ever, our results show that the LS rate only agrees with 
the FGR rate and the full binary collision approach (see 
Ref. in this limit. The differences between the FGR 
and LS results for smaller electron temperatures are due 
to the different definition of "Coulomb logarithms" as 
the FGR approach was shown to reduce to a LS-like for- 
mula 13]. On the other hand, the CM expression (J5]) 
gives considerably reduced energy transfer rates even in 
this high electron temperature limit as the ion accoustic 
modes are well-pronounced here. CM effects only cease 
to exist for high ion temperatures where these modes are 
not present. Thus CM mode effects can be also present 
when describing very hot fusion plasmas. 

The relation between different approaches is studied 
in more detail in Fig. [5l The first panel demonstrates 
how an increasing temperature difference increases the 
lowering of the energy transfer in the CM compared to 
the FGR descriptions. Whereas CM effects slowly de- 
velop with increasing electron temperatures for the lower 
plasmas densities, they are already present at T e =Ti for 
densities above n = 10 23 cm -3 . Here, the electrons are 
degenerate with a Fermi temperature that is higher than 
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FIG. 5: (Color online) Ratio of the CM and FGR energy 
transfer rates for fully ionized hydrogen as a function of the 
electron temperature. Panel a) considers various densities 
while the ion temperature is constant; panel b) shows results 
for different ion temperatures at constant density. All results 
are obtained within the random phase approximation. The 
vertical lines mark the condition T e = Ti. 



the ion temperature and, therefore, ion acoustic modes 
exist even for similar electron and ion temperatures (see 
Eq. © for the occurance of ion acoustic modes). Under 
these conditions, the actual temperature difference be- 
comes irrelevant when searching for CM effects. An ex- 
treme case is given by a density of n — 10 26 cm -3 where 
a lowering of CM rates of up to 90% can be found near 
the maximum of the energy transfer. 

Of course, an increasing ion temperature has a similar 
effect as a decreasing electron density as shown in the 
second panel of Fig. [5] where the well established CM ef- 
fect for Ti = lQ^K vanishes for higher ion temperatures. 
Nevertheless, sufficiently high electron temperatures al- 
ways guarantee that the CM rates are lower than FGR 
or LS rates in these parameter regions. 

Interestingly, the ratio of the CM and FGR rates seems 
to always approach 1/2 in the limit of high electron 
temperatures. In laser-produced plasmas with T e ^$>Ti, 
one should therefore almost always expect reduced en- 
ergy transfer rates due to coupled collective modes. On 
the other hand, the density in a shock-produced plasma 
(T e <Ti) must be high enough to observe CM effects. 

Let us analyze the above mentioned high temperature, 
low density limit of the CM energy transfer rate in more 
detail (see Fig. [5]). The CM and FGR approaches do 
agree in the weak coupling limit if and only if the tem- 
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FIG. 6: (Color online) Ratio of energy transfer rates in CM, 
FGR, and LS in the high temperature (low coupling) limit for 
fully ionized hydrogen. Curves for different ratios of electron 
and ion temperature are given. The three curves below unity 
belong to the ratio of CM to FGR rate. The ratio of FGR 
rate and LS rate is independent of the temperature ratio and 
is given by the green dotted line. The region were the ion 
coupling is small is indicated for each case. 



perature difference between the subsystems is sufficiently 
low so that no coupled collective modes can be excited. 
Larger deviations occur for increased temperature differ- 
ences and smaller electron temperatures (more strongly 
coupled ions; although RPA is used in Fig. [6]). However, 
even for very large electron temperatures the CM effects 
still reduce the coupling between the electron and ion 
components. 

Moreover, FGR and LS have a ratio independent of 
the temperature difference. In addition, numerics show, 
that even in the weakly coupled case with small tempera- 
ture difference, FGR and LS curves intersect rather than 
converge. This is due to the poor cut off in the Coulomb 
integral in LS. Thus the FGR approach rather than LS 
formula should be used to compare with experiments or 
simulations if one searches for coupled mode effects in 
the data. 

Of particular interest is the question whether the easier 
FGR is strictly valid in some cases or not. The difference 
between FGR and CM originates from the different de- 
nominators. Within RPA, one easily obtains 



l (kiu)\ 2 — \e ee {k,uj)\ 2 |e ri (fc,cj)| ; 



+2V ee (k)V u (k) Re xle (*, w) He X % (k, u) 
+2V ee (k)V u (k) Imx° e (fc,w)Im X ? i (fc,w) 



,(k,Lo)f 



1 



1 



(14) 



For the FGR approach to be accurate, the first line only 
most hold. Strictly, this is only possible for very large 
wave numbers k. For the interesting small k values, we 
must require low plasma densities and the static case with 
vanishing imaginary contributions. However, it is essen- 
tial for the occurrence of coupled collective modes and 
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FIG. 7: (Color online) Effect of the (nonequilibrium) ion 
charge state on the energy transfer rates in different approxi- 
mations for a high-Z plasma. Panel a) displays the rates and 
panel b) shows the ratio of the CM and FGR results. The ion 
temperature is Tj = 10 4 K. 



the resulting mechanism of energy transfer to retain the 
frequency dependence in the dielectric functions. When- 
ever, the FGR and CM modes agree (numerics shows 
us that this requires the non-occurance of ion acoustic 
modes), it results from redistributing weights and is a 
numerical and approximate agreement of the expressions. 



B. Charge State and Strong Coupling Effects 

Coupled mode effects on the energy transfer rates are 
of course much more pronounced for plasma ions with 
higher charges since here the ion acoustic modes occur 
for smaller temperature differences or even T e = Tj (see 
Eq. @. These modes exist also for higher k values if 
the ion charge is higher. Different masses of the ions also 
shift the modes to different locations in k space which 
is, however, irrelevant for the (relative) CM effect as the 
modes in FGR behave similarly. 

In addition to the charge state, strong coupling effects 
must be investigated as the CM reduction of the rates 
is stronger for colder system as shown in Fig. [6l In this 
case, the ions can be strongly coupled and not described 
within RPA. Thus, both strong coupling and collective 
excitations must be considered in many situations. 

Figure [7J shows an example for plasmas with highly 
charged ions. When normalized to the electron density 
the energy transfer rates rise linearly with the charge Z in 




T e [K] 

FIG. 8: (Color online) Energy transfer rate for shocked silicon 
at n e = 5.36 x 10 23 cm" 3 , Ti = 54539 K, and with Z = 4 in 
different approximations. The electron Fermi temperature is 
T F = 2.8 x 10 5 K. These conditions are as in the experiment 
of Ng et al. 0. 

the FGR and LS approaches. CM energy transfer rates 
naturally depend stronger on the coupled mode struc- 
ture. Panel b) of Fig. [7J shows that the relative reduction 
for Z = 4 is higher than for Z = 1 so that the charge de- 
pendence of CM rates is smaller than Z. We find charge 
state effects especially for smaller electron temperatures 
where the occurance of the ion acoustic mode strongly 
depends on the charge state of the ion (see in Eqs. ©)■ 
If the electron temperature is already high enough to al- 
low for the mode at Z — 1 , the charge state effect is much 
weaker. For the highest charge state considered in Fig. [JJ 
the influence of the CM rates is even more strongly pro- 
nounced due to strong coupling effects. 

The second example in Fig. [5] considers the situation 
in shocked silicon plasmas as prepared in the experiment 
by Ng et al. @. For this case, the electron tempera- 
ture is lower than the ion temperature. The electrons 
are highly degenerate so that according to Eq. © CM 
effects are possible since the Fermi temperature consid- 
erably exceeds the ion temperature. The ions are also 
strongly coupled, Fa — 64, and local field corrections to 
the response functions have a larger effect as in the ex- 
amples given so far. 

When analysing the case, the first point to realise is 
that neither the usual Landau-Spitzer approach nor the 
Brysk formula [4l| that (partially) accounts for degen- 
eracy effects is well defined. Both formulas contain a 
Coulomb logarithm that either becomes negative or very 
small when considering straight line trajectories or hyper- 
bolic orbits, respectively. Thus, the Coulomb logarithm 
is usually clamped at some larger value. Here, we use 
Ac = 2 for comparison. Still both the LS and Brysk for- 
mula give energy transfer rates that are orders of magni- 
tude larger than the FGR or CM results which acount for 
degenerate electrons in a regular fashion. This is qualita- 
tively consistent with the experimentally observed facts. 
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Comparing the rates obtained by the CM and FGR ex- 
pressions, one makes an interesting observation: on the 
RPA level both approaches yield more or less the same 
results although ion acoustic modes are present in the 
system. The weighting due to the electron Fermi distri- 
butions favors however larger k values where no collective 
mode exist. However, we find a further reduction of the 
CM rates when the effect of strong coupling is included 
by static local field corrections. Since the FGR approach 
naturally includes strong coupling effect [H[ , no changes 
occur within this approach. Accordingly, we find now 
a factor of roughly two between the rates calculated by 
the CM and FGR expressions, respectively. Thus, we 
can confirm the occurance of coupled mode effects under 
these conditions although the main effect arises from the 
correct treatment of degeneracy. 

V. CONCLUSIONS 

We investigated the electron-ion energy transfer rates 
in two-temperature plasmas with special emphasis on 
coupled mode effects. For weakly coupled plasmas, the 
starting point was the quantum version of the Lennard- 
Balescu kinetic equation which naturally includes fully 
coupled collective excitations in the screening function. 
It was shown that a straightforward derivation without 
further approximations yields a coupled mode expression 
for the energy relaxation rates. This kinetic approach 
is consistent and accurate for weakly coupled plasmas 
of any degeneracy. It can be easily extended to plasmas 
with strongly coupled ions by including (static) local field 
corrections. 

For certain conditions, the coupled mode expression 
yield considerably reduced electron-ion energy transfers 
when compared to the LS and FGR approaches. A de- 
tailed analysis showed that this reduction can be traced 
back to the occurance of ion acoustic modes and a re- 
lated redistribution of weight in the dielectric response 
function. This reflects the fact that the coupled mode 
approach considcres a much more realistic mode struc- 



ture. Precise conditions for the occurrence of ion accous- 
tic modes and the coupled mode reduction of the energy 
transfer rates were derived. Interestingly, the CM re- 
duction is preserved for very high electron temperatures 
where the rates are roughly a factor of two less than the 
ones calculated with the FGR or LS approaches. An 
agreement between CM and FGR expressions can only 
be reached if the temperature difference is small or the 
ion temperature is high enough to prohibit the occurance 
of ion acoustic modes. 

The results presented here are concistent with exper- 
imental [1, H, @1 an d theoretical [2(| HH findings of low- 
ered electron-ion energy transfers and thus longer equi- 
libration times. In particular, the strong reduction of 
the equilibration rates in shocked silicon is in qualita- 
tive agreement with measurements. The experiments 
probed however regimes where collective modes as well 
as strong binary collisions and changing potential ener- 
gies are of importance. Thus, a full description must 
also include these effects at least on a level as described 
in Refs. [l(| EH HH . The conditions for the occurance 
of CM effects derived also explain why no or very small 
CM reductions have been seen in recent numerical simu- 
lations mmm. 
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APPENDIX A: DERIVATION COUPLED MODE 
EXPRESSION IN RANDOM PHASE 
APPROXIMATION 

We start with the quantum Lennard-Balescu equation 
for the electron distribution f e (p, t) in homogeneous and 
isotropic plasmas [27| 



dt 



dp' dp dp' 

{2Trhf (2tt?i) 3 {2irh) 3 



Vei(p-p) 



eR pA (p-p,E e (p)-E e (j>),t) 
x2nS (E e (p) + Ei(p') - E e {p) - ^(p')) {^Tif s(p + p' p p') 

X{fe(p,t)fi(p',t) [1 - f e ( P ,t)] [1 - M P ',t)} - f e (p,t)Mp',t) [1 - f e (p,t)] [1 - fi(p' ,t)]} 



(Al) 



r 



The sum runs over all species, but only the ions con- 
tribute to the energy relaxation. E a (p) = p 2 /2m a de- 
notes the kinetic energy of particles of species a and 



Vab(k) — 4:TrZ a Z},e 2 /k 2 the (unscreened) Coulomb poten- 
tial. Dynamic screening is described by the retarded di- 
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electric function which is used here in RPA 

s RPA (p,E,t) = V aa (p) X ° aa (P,E,t). (A2) 

a 

The upper expression relates the screening function to 
the density response function of free particles given by 

XaaKP, ,) J (2ttTi) 3 E + E a (p') - E a (p' + p) + ie' 

(A3) 

which is, in turn, determined by the electron or the ion 
distribution function. 

The electron-ion energy transfer rates are now ob- 
tained by multiplying the Lenard-Balescu equation (jAlj) 
by the electron energy, E e = p 2 /2m e , and an integration 
over the free momentum (see balance equation ([I}. 

For the further proceedings, it is useful to use the 
momentum and the energy transfer during the collision, 
k=p p and Lu — E e (p)—E e (p), respectively. With these 
new variables, one can now apply the relations between 



Fermi-distributions of the plasma species and Bose func- 
tions n B (u>) = [exp(huj/kBT a ) — 1] , namely 

fa(p) [l-f a (p + k)} = 

[Up + k)- Up)] n a B [E a {p) - E a (p + fc)) ,(A4) 

fa(p+k) [l-/ (p)] = 

[fa(p) - /«(p + k)]n B (E a {p + fc) - E a (p)\ . (A5) 



The set of distributions in the third line of the collision 
integral of the Lenard-Balescu equation (lAlj) can then be 
written as 

{/()} = [fe(p + k) - f e (p)] [Up') - Up' + fc)] 

x n B (-uj) n B {uj) - n e B (oj) n B (-ui) . (A6) 

With these transformations, we obtain for the electron- 
ion energy transfer rate 



dp dk doj „ , , , 
E e (p + k) 



«b(-w) n % B (u) - n B {u)n B {-uj) 



V el (fc) 



e R (kuj) 
dp' 



{2nh) 



Up + fc) - Up)] ^s(u -E e {p + k) + E e (p) 
3 [Up') - Up' + fc)] - E^p' + k) + E^p' 



The integral in the second line of this expression is just 
the definition of the imaginary part of the free parti- 
cle density response function for the ionic subsystem 
lmXii{kuj) as defined by (|A3|) . The difference of Bose 
functions in the second line determines the direction of 
the energy transfer. After another variable transforma- 
tion to u>' — —ui, fc' = — fc, and p' =p—k! (changing the 
sign) in the second term proportional to n B (uS) n l B {— u>), 
this second term has the same form as the first one, ex- 
cept that the energy in front of the screened potential is 
E e (p) instead of E e (p— fc). Applying the energy conserv- 
ing (5-function in the collision integral yields 



E e (p + fc) - E e (p) = u. 



(A8) 



Now the remaining electron distributions together with 
the energy conserving (5-function give also the definition 
of an imaginary part of the free density response function; 
this time, the one for the electrons, i.e. Im%° e (S;w). The 
energy transfer rate is thus given by 



0_ 

dt 



E. 



CM 



oc 

= -^E/t5f/ 



duj 
2^ 



V el (fc) 



T RPA 



(ktu) 



xImXe e (M I™-Xii(ku>)n e B (-iv)n B (v) 



(A9) 



(A7) 



r 



If we use the fact that [1/2 — UbU)] and Imx^ a (tj) are 
odd functions with respect to w, we can rearrange the 
upper expression in the form 



oo 

dt (2wh)*j2n 



V ei (fc) 



-RPA 



{kuj) 



x Imx° e (M ImxUkuj) n e B (w) - n l B (uj) 



(A10) 



This expression gives the electron-ion energy transfer rate 
including the effects of fully coupled collective modes in 
the electron-ion system. Similar to the Lenard-Balescu 
equation that was used as a starting point, it describes 
weakly coupled plasmas without restriction with respect 
to the degeneracy. 
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APPENDIX B: EQUIVALENCE TO COUPLED 
MODE DESCRIPTION BY 
DHARMA-WARDANA & PERROT 

The coupled mode expression derived by Dharma- 
wardana & Perrot, equation (50) in Ref. [20l |. reads 



E. 



Tlx 



dk 



du> 



2tt J (2nhy 



AiV CT HA'(fc^)A e (fc,a;) 
|1 Xee(k,w)x«(*,w)l a ' 1 ) 

where we assumed all potentials to be of Coulomb type 
for simplicity. The difference of Bosc functions is here 
the same as in Eq. (|A10[) : AN ei = n B (to) — n B {u). The 
functions A a (k,uj) are related to the imaginary parts of 
the response functions Xaa(k,(jj) 



A a {k,u) = -2lm X aa(k,uj). 



(B2) 



One should however notice that these are full density 
response functions of a coupled system. 

The density response function Xaa for a single species 
can also be written in terms of the polarization function 
II aa and the dielectric function e of the (fully coupled) 
medium. Neglecting cross terms, we have [421 ] 



Xaa(k 7 U)) 



_ Tlaa{k,Uj) 

1 - V aa (k)Tl aa (k, u) e aa (k,u) 



U aa (k,uj) 



(B3) 



This is nothing else than the equation of motion for the 
density response function. Real and imaginary parts of 
the response function can then be expressed as 

n (l . \ Re£q(fc,a;) - \e a {k,uj)\ 2 
Rexaa(fc,w) = ~ ~ j-^ n _ ^ , (B4) 



lmxaa(k,uj) = - 



V aa (k)\s a (k,u)\ 2 
Ime a (k,uj) 

14a(fc)|£a(fc,^)| 2 " 



(B5) 



With the help of these expressions, the denominator in 
expression (|B1I) which constitudes the differences to the 
FGR formula becomes 



|l-V£(fc)Xee(fc,w)X«(fc,w)| S 



\e(k,u)f 



\^e(k,Uj)\ 2 \Ei(k,L0)\ : ' 



(B6) 



The dielectric function in the numerator of Eq. (|B6[) is 
the one for the full system given by e = 1 — ^ Q V aa Il aa . 
The partial dielectric functions e a in this expression are 
cancled by the one contained in the imaginary part of 
the full density response function Xaa when the latter 
are written as in eq. (|B5|) . Accordingly, one obtains an 
expression identical to Eq. (|A10p when inserting these 
definitions and transformations into Eq. (|Bip . 



APPENDIX C: CONDITIONS FOR THE 
OCCURANCE OF ION ACOUSTIC MODES 

It is well known that ion accoustic modes occur in non- 
equilibrium plasmas with T e ^> Ti [37|]. Here, we give a 
more precise estimate of the necessary temperature rela- 
tion needed to allow for this type of collective modes in 
weakly coupled plasmas. 

For ion acoustic modes to exist, the real part of the 
total dielectric function must vanish. We first consider 
the ionic contribution in non-degenerate limit which is 
almost always sufficient for the heavy ions. In RPA, the 
polarisation function is identical to the free particle den- 
sity response function, that is = Xu- ^ s rea ^ P ar t 
can be written in terms of the confluent hypergeometric 
function [13, 51 



Be X %(k,iv) = 



knT 



1 



p 2 k B T 



iFi 



4 



2p 2 k B T 

( C1 ) 

which is defined over an integral [431 ]. However, one can 
approximate this function for the special first two argu- 
ments needed by the following Pade formula (2?| 



(C2) 





"*) 










14 


-§ + 


- — ^ 

10 n 


42 ~ 


218 ~ 


7x 5 +x e 
9360 




X 2 

1 


r 6 


+ Xi - 
' 24 


1 120 


, X 6 , X 7 

1 720 + 4860 



As a required condition for ion acoustic modes to oc- 
cur, the minimum of the response functions (|Clj) must at 
least compensate the electronic contribution plus unity. 
We thus search for an approximation of Eq. (jC2[) that 
conserves this minimum. Performing a frequency deriva- 
tive, we find its location to be approximately at u>o = 
2.36 k (ksTi/mi) 1 ^ 2 . Using the minimum frequency ujq 
in Eq. (|C1[) with the hypergeometric function given by 
Eq. IJC2J) yields 



Ree»(fc,wo) = 1 - 0.27 



k 2 



(C3) 



where Ki — {AirZ 2 e 2 ni/kBTi) 1 l 2 is the inverse of the ion 
part of the classical Debye screening length. 

The electronic part of the dielectric function can be 
very well estimated in the zero frequency (static) long 
wave length limit 



Ree e (fc,0) = 1 



k 2 



(C4) 



Again, k denotes the inverse screening length. For small 
momenta, the unity in Eqs. (|C3j) and (|C4[) can be ne- 
glected. The condition for the ion acoustic modes to 
emerge becomes then 



k 2 < 0.27 k? 



(C5) 



E 



(|C5[) quantifies the long known condition of Tj <C T e 
for the ion acoustic waves to occur. Direct relations 
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for the electron and ion temperatures needed may be 
obtained by inserting an appropriate expression for the 
electron screening length. In the nondegenerate limit, 
it is given by the electron Debye length and one ob- 
tains the upper relation of Eq. ([9]) . For highly degener- 
ate electrons, Thomas-Fermi screening with n e = k^f = 
4m e e 2 (3n e /7r) 1 / 3 must be used. Since the Thomas-Fermi 



length is temperature independent, the role of the elec- 
tron temperature is taken over by the electron Fermi tem- 
perature which yields the second condition in Eq. ([9|) 
for the existence of the ion acoustic modes. Accord- 
ingly, they can also occur in degenerate plasmas at equi- 
librium or even in shock-produced high-density plasmas 
with T e <Ti. 
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